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ABSTRACT 

The Space Infrared Telescope Facility (SIRTF) will revolutionize the study of dust-obscured star for- 
mation in distant galaxies. Although deep images from the Multiband Imaging Photometer for SIRTF 
(MIPS) will provide coverage at 24, 70, and 160/im, the bulk of MLPS-detected objects may only have 
accurate photometry in the shorter wavelength bands due to the confusion noise. Therefore, we have ex- 
plored the potential for constraining the total infrared (IR) fluxes of distant galaxies with solely the 24/im 
flux density, and for the combination of 24/im and 70/im data. We also discuss the inherent systematic 
uncertainties in making these transitions. Under the assumption that distant star-forming galaxies have 
IR spectral energy distributions (SEDs) that are represented somewhere in the local Universe, the 24/im 
data (plus optical and X-ray data to allow redshift estimation and AGN rejection) constrains the total 
IR luminosity to within a factor of 2.5 for galaxies with 0.4 < z < 1.6. Incorporating the 70/im data 
substantially improves this constraint by a factor < 6. Lastly, we argue that if the shape of the IR SED 
is known (or well constrained; e.g., because of high IR luminosity, or low ultraviolet/IR flux ratio), then 
the IR luminosity can be estimated with more certainty. 

Subject headings: cosmology: observations — galaxies: high-redshift — infrared: galaxies 



1. INTRODUCTION 

The evolution of the global, volume-averaged star for- 
mation rate (SFR) is currently a topic of intense interest 
(see, e.g., Madau et al. 1996; Steidel et al. 1999; Yan et 
al. 1999; Haarsma et al. 2000). Observationally, a galaxy's 
SFR must be inferred from the luminosity density at wave- 
lengths dominated by young stars, such as the ultraviolet 
(UV), nebular emission lines, mid-to-far infrared (IR), or 
radio (see, e.g., Kennicutt 1998; Condon 1992), and by 
making assumptions about the form of the IMF. Because 
of dust, rest-frame UV/optical indicators often underes- 
timate the intrinsic SFR, especially for more luminous 
galaxies (see, e.g., Calzetti, Kinney, & Storchi-Bergmann 
1994; Calzetti 2001; Bell 2003). Based on the UV/optical 
spectral energy distributions (SEDs) and far-IR emission 
of distant galaxies and QSOs, dust appears to be prevalent 
in high-redshift objects (z ~ 2 — 5, e.g., Pettini et al. 1998; 
Adelberger & Steidel 2000; Carilli et al. 2000; Papovich, 
Dickinson, & Ferguson 2001; Chapman et al. 2002). Dust 
reprocesses the energy absorbed from the UV/optical into 
the mid/far-IR, and thus galaxy IR luminosities allow one 
to 'balance the energy budget' and paint a much more com- 
plete picture of star formation throughout cosmic history 
(see, e.g., Sanders & Mirabel 1996; Blain et al. 1999; Flores 
et al. 1999). 

The Space Infrared Telescope Facility (SIRTF), to be 
launched in 2003, offers to revolutionize our view at IR 
wavelengths with unprecedented sensitivity and resolu- 
tion. Several cosmological surveys are being planned to ob- 
tain 24, 70, and 160/Ltm data with the Multiband Imaging 
Photometer for SIRTF (MIPS; see, e.g., Lonsdale 2001; 
Dickinson & Giavalisco 2001; Rieke 2001; Rieke et al. 
2001). The combination of these three bands will span the 
peak of the IR SEDs and provide a fairly robust tracer of 

1 Present Address: Max Planck Institut fur Astronomie, Konigstuhl 
17, D-69117 Heidelberg, Germany; bell@mpia.de 



the total IR emission (e.g., Dale & Helou 2002). However, 
MIPS observations will rapidly become confusion limited 
due to the increasing source density at faint fluxes (e.g., Xu 
et al. 2001; Dole, Lagache, & Puget 2002), especially for 
the longer wavelength data, which have lower sensitivity 
and resolution (MIPS resolution is roughly proportional to 
the bandpass central wavelength). Thus, full coverage of 
galaxies' IR SEDs with accurate photometry may only be 
possible for relatively nearby or bright objects. 

In this Letter we consider the uncertainties inherent in 
translating MIPS photometry into total IR fluxes in the 
case that object photometry is only available from the 
shorter wavelength bands. The relation between mid-IR 
and far-IR is complex. Thus we explore the connection 
between the 24/xm data and the total IR emission from 
both the observational (§2) and modeling (§3) perspec- 
tive, with the assumption that local galaxy SEDs are rep- 
resentative of high-redshift analogs (which has not been 
robustly demonstrated, although see Adelberger & Steidel 
2000; Elbaz et al. 2002). In §4, we investigate improving 
these constraints using 70/im data, the galaxy luminosity, 
and/or UV flux. We summarize our results in §5. Because 
extensive X-ray /optical/near- IR coverage will be available 
for the majority of deep MIPS survey regions (see refer- 
ences above) , we assume here that galaxies with dominant 
AGN contributions can be rejected using the X-ray data, 
and that galaxies' redshifts (spectroscopic or photometric) 
will be known. 

2. OBSERVATIONAL CONSTRAINTS BETWEEN THE MID-IR 
AND TOTAL IR EMISSION: THE CASE AT Z ~ 1 

To explore the correlation between mid-IR fluxes and 
the total IR flux 2 for local galaxies, we have taken data 
from the Infrared Astronomical Satellite (IRAS) at 12/im. 

2 We define the total IR flux as the integrated emission from 
8 — 1000/im, i.e., Fm = F(S — 1000/im), estimated using direct in- 
tegration from 8 — lOO^m and extrapolated to 1000/im using a A -1 
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Fig. 1.— The ratio of the IRAS 12//m to total IR flux (from 
8 — 1000pm) as a function IRAS 100pm to 60pm flux-density ra- 
tio; see Bell (2003) for details. Symbols denote various galaxy types: 
starbursting galaxies (open stars); ultra— luminous IR galaxies (Bllcd 
stars); normal galaxies (open circles). The error bar indicates typ- 
ical uncertainties. Results for model galaxies arc shown as lines 
(with line definitions as in figure 2). 

Note that the MIPS 24^im band in the observed frame for 
z = 1 nearly matches the rest-frame IRAS 12 rim band, 
with a if-correction that varies weakly with galaxy SED: 
Iv(MIPS24pm, z = 1) « 1.3 L V (IRAS 12fim, z = 0). 

Figure 1 shows the IRAS 12pm to total IR flux ratio as 
a function of IRAS 100pm to 60pm color for 156 galaxies 
spanning all star-forming types (Bell 2003). 3 The figure 
illustrates that a nearly constant proportion of the total 
IR flux emerges in the mid-IR for local galaxies. The mean 
ratio is F(12pm)/i<iR = 0.16, with a la scatter between 
0.10 and 0.22. This is consistent with IRAS and ISO re- 
sults, which show a relatively constant conversion between 
rest-frame mid-IR (7/im— 15/xm) and total IR luminosities 
(e.g., Spingolio et al. 1995; Chary & Elbaz 2001; Roussel 
et al. 2001). Clearly, in local galaxies the dust responsible 
for the unidentified IR bands (UIBs; most plausibly associ- 
ated with polycyclic aromatic hydrocarbons (PAHs), e.g., 
Puget & Leger 1989; Genzel & Cesarsky 2000) that domi- 
nates the flux in the mid-IR correlates fairly well with the 
dust that produces the bulk of the total IR flux emission. 

Importantly, figure 1 shows that the ratio of 12pm to to- 
tal IR flux does not sensitively depend on the temperature 
of the dust that dominates at far-IR wavelengths (as mea- 
sured by /„(100pm)//„(60pm)): the systematic variation 
with dust temperature seems to be smaller than the scat- 
ter. Overlaid in the figure are the results of models from 
the literature with different treatments of the mid/far- IR 

cmissivity. These fluxes are typically a factor of 2 larger than those of 
Hclou et al. (1988) and 30% higher than those of Sanders & Mirabel 
(1996), and are consistent to within 10% with values extrapolated 
using 8-170pm data from IRAS and Tuffs et al. (2002). 

3 We denote the average flux density measured from a bandpass 
with central wavelength A as /^(A), and we define the flux from the 
bandpass as F(X) = v f v (\), where v is the frequency corresponding 
to the wavelength A. 



SEDs and that span a range of dust temperatures and IR 
luminosities (described in §3). These models broadly re- 
produce the mean and scatter observed in the plot. These 
IRAS data suggest that if the relationship between the 
mid-IR and total IR emission at high-redshifts is similar 
to that observed locally, then the MIPS 24pm data will 
provide an estimate to the total IR emission for z <~ 1 
galaxies. Clearly, the scatter is significant (a factor of 
« 2.5), and is largely systematic in the sense that it can 
only be improved upon by better constraining the shape 
of the galaxy IR SED. 

3. MODELING THE MID-IR TO TOTAL IR EMISSION: 
EXTENDING THE REDSHIFT RANGE 

Here, we turn to models to test the result from §2 for 
a wider range of rcdshifts. In order to gauge the system- 
atic model uncertainties, we explore three independently 
constructed sets of models, i) We use the STARDUST tem- 
plates for 17 local galaxies with a range of IR luminosities, 

5 x 10 9 < Lir/L q < 4 x 10 12 , from Devriendt, Guiderdoni, 

6 Sadat (1999). Note that the luminosities of the seven 
Virgo Cluster galaxies in this model were underestimated 
by a factor of between 10 and 300 (J. Devriendt, 2002, 
private communication) ; this does not affect our results as 
we use only colors (flux ratios), it) We have also used the 
fits of the grasil model of Silva et al. (1998) to six galax- 
ies (Arp 220, M51, M82, M100, NGC 6090, NGC 6946). 
Hi) Lastly, we use the empirical models based on local 
galaxy ISO data from Dale et al. (2001) with a range 
of far-IR flux ratios (a surrogate for dust temperature), 
-0.4 < logF(60pm)/F(100pm) < 0.5. These models ap- 
ply explicitly to local galaxy SEDs, and thus the caveat 
persists that they may not describe high-z galaxies. 

To derive the observed-frame mid-IR flux as a func- 
tion of redshift, we averaged the model galaxy SEDs with 
the MIPS 24pm bandpass (including effects from the fil- 
ter transmission and detector response) over a range of 
rcdshifts. In figure 2, we show the ratio of the observed- 
frame MIPS 24/im to rest-frame total IR flux versus red- 
shift. As expected from the initial analysis of 12pm IRAS 
data, the MIPS 24pm emission is a fair tracer of the total 
(rest-frame) IR emission, especially for the redshift range 
0-4 S z iS 1-6, for which the scatter is a factor of 2 — 3. 

It is worth discussing this in somewhat more detail. At 
z <~ 0, the observed 24pm band traces emission from UIBs, 
and very small grains that have been transiently raised 
to high temperatures by absorption of single UV photons 
(Desert, Boulanger, & Puget 1990). Model SEDs show 
great variation in the prominence of this component, de- 
pending largely on the relative intensity of the UV ra- 
diation field in a given galaxy. As one moves toward 
z ~ 1 — 2.5, the MIPS 24/im data probes lower rest-frame 
wavelengths where the UIBs increasingly contribute to the 
flux especially in discrete features from ~ 7 — 13pm (Gen- 
zel & Cesarsky 2000, and references therein). Empirically, 
these UIBs are rather more prominent for 'cooler' galaxies 
with less intense UV radiation fields (causing the tracks for 
warmer and cooler galaxies to cross at z ~ 0.5). The 'dip' 
in 24pm-to-total IR ratio at z ~ 1.4 is attributed to gaps 
between the UIBs. Silva et al. (1998), Devriendt et al. 
(1999), and Dale et al. (2001) model the UIBs somewhat 
differently. The largest difference between the models at 
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Fig. 2. — The ratio of the observed MIPS 24 /im emission to 
total (rest-frame) IR emission as a function of redshift. Curves 
display the models described in the text (labeled in the figure inset). 
The arrow depicts galaxy SEDs from "Warmer" to "Colder" dust 
temperatures. 

z > 1 is the contribution of the stellar continuum to the 
mid-IR emission. The galaxy templates match well for 
warm galaxies, and diverge for cooler galaxies, as Dale et 
al. do not have a large stellar contribution in their cold 
galaxies, whereas Devriendt et al. include seven Virgo 
Cluster galaxy templates that are dominated by stellar 
emission at A < 15/im (Boselli et al. 1998), and Silva et al. 
study three quiescently star-forming spirals (M100, M51, 
NGC 6946) that contain some stellar contribution to the 
mid-TR emission. Certainly, the contribution of the old 
stellar population to dust heating is correlated with dust 
temperature (see, e.g, Bell 2003, and references therein); 
however, the scatter is large, leaving the issue of stellar 
contamination of A < 15/zm fluxes largely open. More- 
over, how these processes contribute to the IR SEDs of 
high-z galaxies is almost entirely unknown. 

Given this (and with the inherent assumption that dis- 
tant star-forming galaxies have IR SEDs that are repre- 
sented somewhere in the local Universe), one expects that 
at 0.4 < z < 1.6 it should be possible to constrain the 
total IR flux to within a factor of 2 — 3 using 24/im data 
alone: i 7 '(24/im) represents 5 — 25% of the total IR flux. 
Outside this redshift range, the conversion is rather less 
certain. A factor of « 2 — 3 in error in total IR flux may be 
quite adequate, but can potentially frustrate many scien- 
tific analyses. Because this uncertainty is largely system- 
atic, one must include additional knowledge of the shape 
of the galaxy SED in order to improve the constraint on 
the total IR flux. 

4. IMPROVING THE TOTAL IR EMISSION CONSTRAINTS 

Based on the models illustrated in figure 2, the scat- 
ter between the 24/im flux and the total IR flux emission 
stems primarily from a lack of knowledge on the shape 
of the galaxy IR SED (and/or dust temperature). MIPS 
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FlG. 3.— The ratio of the observed SIRTF /MIPS 24 /im emis- 
sion to the total (rest-frame) IR emission as a function of MIPS 
24/^m/70/^m flux ratio. Curves with symbols illustrate the various 
predictions from the models described in the text averaged over 
the redshift interval (as listed in the figure inset), with line styles 
as indicated in figure 2. The thick solid and dotted lines show a 
second-order polynomial fit and its standard error to the average 
over all model predictions and rcdshifts (see text). 

70/im flux measurements will be available for a signifi- 
cant fraction of MIPS 24/im sources (although the exact 
fraction and redshift distribution are strongly model de- 
pendent), which should enhance the constraints on the 
total IR emission. Here we explore the constraints on 
the total galaxy IR emission in the case that only MIPS 
24/im and 70/im data are available. In figure 3 we show 
the ratio of the observed MIPS 24/im to total IR emis- 
sion versus the ratio of the observed MIPS 24 and 70/im 
fluxes in different redshift intervals. There is generally 
very little scatter in the relationship for these redshifts 
and colors: the ratio of the 24-to-70/im fluxes seems to 
correlate with dust temperature over a wide range of red- 
shifts. To crudely parameterize the relationship between 
the IR fluxes in these bandpasses, we have fit a second- 
order polynomial to the average over all models for all 
redshifts shown in figure 3. We find a best fit to the equa- 
tion, F(24/*m)/F(8 - 1000/im) = A + Bx + Cx 2 , where 
x = logF(24/im)/i 7 '(70/im), with polynomial coefficients, 
A = 0.160, B = 0.135, and C = 0.033, and a standard 
error of <5[F(24/im)/-F(8 - 1000/im)] = 0.036. This fit and 
error are overlaid in figure 3. Our parameterization pro- 
vides a means for quantifying the scatter in estimating the 
total IR flux emission using the MIPS 24 and 70/im data 
(although the physical relation between these fluxes and 
the total IR emission is of course more complicated than 
this simple formalism), and it spans the systematic uncer- 
tainties inherent in the model predictions and is largely 
independent of redshift between 0.3 < z < 3.3. Therefore 
with both MIPS 24 and 70/im detections (and with the as- 
sumption that local galaxy IR SEDs apply to high-redshift 
analogs) this result suggests that the uncertainty on the 
total IR flux is improved by up to a factor ~ 6 (relative to 
that achieved using 24/im data only). 
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Lastly, if distant luminous IR galaxies (LIRGs; Lm > 
lO n L0) have comparable SEDs to those found locally, 
then it should be possible to better constrain the IR lu- 
minosity for this class of galaxy. Locally, LIRGs have 
"Warmer" dust temperatures (Sanders & Mirabel 1996), 
thus allowing for a more robust conversion between ob- 
served 24um and total IR luminosity (adopting the warmer 
tracks, with less stellar contribution, in figure 2). It may 
be possible to classify these galaxies as LIRGs on the ba- 
sis of their high IR luminosities (indeed, although model 
dependent, it may be the case that only LIRGs will be 
detectable in confusion-limited surveys for z > 1; e.g., 
Dole et al. 2002). Alternatively, many local LIRGs have 
low UV-to-IR ratios F(1550A)/fi R < 0.02, which can be 
derived from observed-frame [/-band to 24/itm ratios for 
z ~ 1 galaxies, albeit in a model-dependent way. These 
conversions obviously depend sensitively on high rcdshift 
LIRGs acting in similar ways to local LIRGs, which of 
course must be tested using spectral observations along 
with 70 and 160/mi photometry for the brightest galaxies. 

5. CONCLUSIONS 

SIRTF will greatly improve our understanding of the 
IR emission of (and hence the star-formation processes 
within) distant galaxies. Because deep surveys could un- 
cover many star-forming galaxies in the shorter-wave- 
length MIPS data with no longer-wavelength counter- 
parts, we have explored the efficacy of constraining the 
total IR galaxy emission using the data at 24/im only and 
the combination of 24 and 70/im (and under the assump- 
tion that ancillary optical and X-ray data are available to 
allow redshift estimation and AGN rejection). Assuming 
that distant star-forming galaxies have IR SEDs that are 
represented somewhere in the local Universe, the 24/im 
data should constrain the integrated IR flux to within a 
factor of 2.5: F (24/j,m) / F m ~ 5 - 25% for galaxies with 
0-4 < z < 1.6. Including MIPS 70/iin data, the IR lumi- 
nosity can be estimated with considerably more certainty 
(a factor < 6 improvement over the uncertainties for using 
24/im data only). Lastly, if one can assume that the shape 
of the galaxy SED is similar to local luminous IR galaxies 
(L > 10 11 Lq\ e.g., because of high IR luminosity or low 
UV/IR flux ratio), then the constraints should improve. 

Throughout this study we have made the explicit as- 
sumption that the relationships between the mid- and 
total-IR emission and the shape of the IR SEDs observed 
locally apply to properties of high-redshift galaxies. As a 
final caveat, we emphasize that for many reasons this may 
not be the case. Indeed, the population of IR-emitting 
galaxies must undergo strong evolution in order to match 
the cosmic IR background (Hauser & Dwck 2001; Chary 
& Elbaz 2001; Xu et al. 2001; Elbaz et al. 2002; Dole et 
al. 2002), and it is unclear how such evolution manifests 
itself. The mid IR UIBs (e.g., PAHs) and total IR emis- 
sion relation at high redshift may be different due to, e.g., a 
changing dust composition, or a significantly different dust 
heating from older stellar populations than that observed 
locally. Moreover, we have assumed that the shape of the 
IR SED as a function of IR luminosity for local galaxies 
applies to high-redshift analogs, which is poorly known 
(although not inconsistent with observations; see, e.g., 
Adclberger & Steidel 2000). Only with large samples of 



galaxies detected in all three MIPS bands, and with spec- 
troscopic measurements from the SIRTF Infrared Spec- 
trograph (IRS), will these assumptions be testable. Once 
such constraints are established, the data should provide a 
prescription for relating the MIPS data to total IR fluxes. 
Given an estimate on the total IR luminosity, it is then 
possible to estimate SFR (modulo the usual sources of un- 
certainty; see, e.g., Kennicutt 1998; Bell 2003). 
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